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Searches for supersymmetric top quarks at the LHC have been making great progress in 
pushing sensitivity out to higher mass, but are famously plagued by gaps in coverage around 
lower-mass regions where the decay phase space is closing off. Within the common stop- 
NLSP / neutralino-LSP simplihed model, the line in the mass plane where there is just 
enough phase space to produce an on-shell top quark remains almost completely uncon¬ 
strained. Here, we show that is possible to dehne searches capable of probing a large patch 
of this difficult region, with S/B ~ 1 and signihcances often well beyond 5a. The basic 
strategy is to leverage the large energy gain of LHC Run 2, leading to a sizable population 
of stop pair events recoiling against a hard jet. The recoil not only re-establishes a 
signature, but also leads to a distinctive anti-correlation between the and the recoil jet 
transverse vectors when the stops decay all-hadronically. Accounting for jet combinatorics, 
backgrounds, and imperfections in measurements, we estimate that Run 2 will already 
start to close the gap in exclusion sensitivity with the hrst few 10s of fb“^. By 300 fb“^, 
exclusion sensitivity may extend from stop masses of 550 GeV on the high side down to 
below 200 GeV on the low side, approaching the “stealth” point at = mt and potentially 
overlapping with limits from tf cross section and spin correlation measurements. 



I. INTRODUCTION 


Light stops with mass below a TeV are extremely well-motivated by the supersymmetric 
solution to the hierarchy problem. The uniquely important role of these particles has inspired 
a growing and increasingly sophisticated set of dedicated searches at the LHC, targeting an 
array of different possible decay topologies HHIZI (see also HSHg). While these searches 
have already probed signihcant portions of the possible model space below a TeV, sizable 
gaps in coverage remain even at 0(100 GeV), leaving us to consider: Is it possible that 
light stops have already been produced in abundance in LHC Run 1 but have simply been 
missed? 

In perhaps the most minimalistic benchmark scenario, stops are produced directly in pairs 
via QCD, and each stop undergoes a one-step R-parity-conserving cascade into an invisible 
neutralino LSP and an on-shell or off-shell top quark: 

pp—^ii*, f —)■ -t-(1 ) 

The visible composition of the hnal state is then identical to that of ti, which serves as 
a copious background. The main kinematic handle exploited in most searches has been 
the additional injection of (or more properly from the neutralinos. For S> m^, 
exclusion limits from ti -|- searches at Run 1 extend beyond 700 GeV [10]. However, 
such searches face a major challenge when confronted with lower-mass regions in the stop- 
neutralino mass plane where the is squeezed out. In particular, much attention has 
recently been directed at the “top compression line” ~ + rrit, which dehnes the 

boundary between two-body decays into an on-shell top quark and neutralino, and three- 
body decays via an off-shell top quark into Wbx^. Limits along this compression line are 
largely nonexistent over a roughly 20 GeV-wide gap in stop mass. 

Proposals to probe this region using the total ti cross section and spin correlations [20l 
\2T\ have led to some inroads near the so-called “stealth” point (m^, m^) = (mj, 0 ) [U 0 ], 
but theoretical limitations make it unclear if these searches can be pushed much further. 
The relatively long lifetimes of stops very near to the top compression line has led to a 
complementary suggestion to use the annihilation-decays of stoponium [ 221 - I 21 ] . which would 
lead to distinctive resonant diboson signatures (including, e.g., 77 and Z^)- Projections for 
Run 2 predict sensitivity up to stop masses of several hundred GeV, depending in detail on 
the stop chirality admixture. However, these searches become insensitive if the individual 
stops decay more quickly than the stoponium, which generally occurs as soon as the stop- 
neutralino mass difference opens up to even O(GeV). Other approaches have sought to use 
the small amount of that is available within the bulk of the produced stop pair events. 
Very detailed measurements of the shapes of the tails of .^^-sensitive observables may be 
promising [25], but a careful accounting of theoretical and experimental errors is not always 
available, and the one measurement of this type that has been carried out [ 2 ] (by ATLAS, 
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in the /+jets channel) does not reach the compression line. A simple cut-and-count style 
search based on dileptonic mT 2 [261 EZ! should still be viable due to a particularly sharp 
turnoff of the background, and is also sensitive near the stealth point. But the maximum 
mass reach of such a search is ultimately limited by low statistics and exhibits a significant 
dependence on stop chirality. It has also been suggested to utilize electroweak production of 
stop pairs via VBF, albeit with difficulties in probing stops much heavier than rrit PSI1^ . 

Given these various limitations, there remains a clear need to consider further alternative 
options, lest comprehensive exploration of the top compression line be deferred to future 
precision lepton colliders. To make progress, we may take some inspiration from another 
compression line, at the very lowest end of the stop mass range: mi ~ m^. There, not 
only the but all visible activity is being squeezed out of the decay. Nonetheless, limits 
exist from the LHC, presently up to roughly 260 GeV mm- These are obtained using the 
classic trick of cutting into the region of production phase space where a sparticle pair is 
produced in association with a visible hard recoil particle, in this case a jet. For an almost 
completely compressed spectrum, the neutralinos go to zero velocity in the rest frames of 
their parents, but carry the full energy and therefore take up the full four-momenta. For 
stop pair production, the vector in lab-frame is then automatically equal to the net it* 
transverse momentum vector, which in turn approximately balances against the leading jet. 

For stop-neutralino spectra near the top compression line, we can define an analogous 
trick, but now face several novelties. The neutralinos again approach zero velocity in 
their parent frames, but they share the four-momenta with (almost) on-shell sister top 
quarks, with fraction m^/mi taken up by the neutralinos. Therefore, in the limit of per¬ 
fectly compressed two-body decay i —)■ and assuming a single dominant recoil jet with 

pV(jet) ~ —pTiii*), we get the following relation, 

- -Pr(jet) X (2) 

mi 

The is now attenuated relative to the recoil pr, by a factor that can nominally extend 
down to zero in the massless neutralino limit (corresponding to the stealth stops [20]). This 
attenuation will generally make searches much more challenging when -C m^ along the 
compression line, such that great care will be required in understanding the lower mass 
reach. For a given neutralino mass, the extra 2mt worth of energy required to make a 
stop pair also leads to much lower rates relative to conventionally compressed spectra with 
mi ~ m^, especially in association with a proportionately energetic recoil jet. This issue in 
particular will be greatly ameliorated with the higher beam energy of the upgraded LHG. 
Finally, the two stop decays produce two on-shell or off-shell top quarks, which add to the 
visible activity and can inject further if either W decays leptonically. Perhaps somewhat 
counterintuitively, the cleanest signal may then be the all-hadronic decay mode, where all 
of the comes from the neutralinos, and Eq. is most closely followed. However, this 
decay mode also maximizes possible QGD backgrounds, as well as our possible confusion 
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over exactly which jets come from the recoil against the stop pair versus from their decays. 

The possible utility of high-momentum recoils in this respect was emphasized relatively 
recently in [30] . In the present paper, we seek to put these ideas on firmer phenomenological 
footing, including a novel set of cuts and treatment of jet combinatorics, a detailed ac¬ 
counting of the various backgrounds, and allowance for a range of possible measurement 
performances. Targeting all-hadronic stop decays, we typically find a healthy S/B ~ 1, 
ensuring robustness against systematic errors of up to 0(10%). We proceed to make a 
detailed forecast for the possible discovery and exclusion coverage in the stop-neutralino 
mass plane. Our results are summarized in Fig. where the proposed search is seen to 
cover a large portion of the formerly inaccessible top compression line, acting as a bridge 
between the two-body and three-body search strategies. For the expected 300 fb“^ to be 
delivered through Run 3 of the LHC, exclusion sensitivity extends up to 550 GeV. On the 
lower end, shrinking poses a major complication, but we find that exclusion sensitivity 
down to ~ -|- 0(10 GeV) may be possible. This would merge our forecasted coverage 

with that of ti cross section measurements and other techniques that perform well in the 
stealth region, allowing for unbroken coverage. If this can be achieved, it would be a major 
accomplishment of the LHG, and a further demonstration that the enormous luminosity 
and broad bandwidth of accessible energies there provides unique opportunities, even for 
relatively low-mass physics with subtle kinematics. 

Our paper is organized as follows. The next section outlines our proposed analysis strat¬ 
egy and presents our estimated signal sensitivities. Section 0 discusses the results and 
possible extensions. More details of the generation of our event samples are presented in 
Appendix 


II. PROPOSED ANALYSIS AND PREDICTED COVERAGE 

Our proposed analysis requires only a few ingredients: 

• A veto on isolated leptons. 

• A high multiplicity of jets and at least two 6-tags. 

• An energetic “ISR-jet” candidate. 

• Coarsely-reconstructed top-candidates whose masses are not significantly above m*. 

• A strong anticorrelation of ISR-jet and directions. 

• A “significant” amount of missing energy, ^t/\/H r, localized near a value set by the 
ISR-jet pt cut and m^/mp 
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In more detail, our full reconstruction and selection, applied to 13 TeV simulated data 
(Appendix]^, proceeds as follows. 

Reconstructed electrons (muons) are first selected starting from truth leptons with pt/ > 
10 GeV and \rii\ < 2.5, and flat identification efficiency of 0.90 (0.95). Electrons are then 
isolated by first computing Yhi bT,i|AJ?<o .2 (where the sum is over all other particles within 
AR < 0.2 of the electron) and requiring 

^ I Ai2<0.2 ^ Q (3) 

PT,e 

Electrons that fail this isolation criterion, as well as all other unidentifiable leptons, are 
returned to the particle list as “hadrons” to be used in jet clustering. Additionally, there 
must be no jets (defined below) within 0.4 of either an electron or muon. Otherwise, the 
lepton is vector-summed into the closest jet.^ Events that contain any surviving isolated 
leptons are then discarded. This lepton veto signihcantly reduces important backgrounds 
where the arises from a W boson decay, especially /-|-jets ti events and leptonic hh-l-jets. 
More aggressive approaches than ours are also possible, using r anti-tagging and/or vetoes 
on more loosely-identified leptons. Ultimately, we find that our backgrounds containing Ws 
are moderately dominated by rz/,-. 

Jets are clustered from all truth hadrons, photons, and unidentihed leptons (including 
electrons that fail the initial isolation step). The anti-fc^ algorithm [21] in Fast Jet [32| 
is applied with R = 0.4, an initial px threshold of 15 GeV, and \p\ < 5.0. Jets 
from this stage are used for the lepton isolation above. Individual jet energies are then 
smeared with gaussians according to the expectation for the Run 2 & 3 conditions of 
~ 50 simultaneous pileup events, as projected in the Snowmass 2013 simulation note [33] : 
<7{pt)/pt = (8.2 GeV)/pT © (0.55 GeV^^^)/y/pT © 0.02.^ Subsequently, an event must have 
at least seven reconstructed jets with smeared pt above 20 GeV, highly favoring the all- 
hadronic tU+jet signal topology and further reducing backgrounds.^ 

Jets with Ir^l < 2.5 are 5-tagged according to an assumed working point with an efficiency 
of 0.70 (0.10) for truth &-jets (c-jets). Jets are first truth ffavor-tagged by looking for the 
heaviest overlapping b- or c-hadron in the event record, and then assigned a reconstruction- 


^ While these steps do not explicitly fold in pileup, significant drops in lepton reconstruction and isolation 
efficiencies in the coming LHC runs are unlikely, especially given the availability of isolation methods that 
are more tracker-based. It is also important to note that, because of the high recoil pr cut demanded 
below, leptons in the dominant backgrounds tend to be quite energetic. 

^ As of this writing, the most recent version (vl) contains a shifted-decimal typo for the noise coefficient in 
the written formula. 

^ We do not model “pileup jets” consisting mostly of diffuse pileup particles, of which 0(2) per event are 
expected [JJ] given our px threshold and before dedicated pileup-jet rejection. We anticipate that these 
will be rejected with reasonable enough efficiency (see, e.g., [35]) so as not to have a major impact on our 
analysis, though higher thresholds on the individual jet pxs would also be an option if necessary. 
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lepton veto 

no isolated ID’ed leptons with pril) > 10 GeV, \r]{l)\ < 2.5 

jets 

PtU) > 20 GeV, \r]{j)\ < 5; N{j) >= 7, A(6-tag) >= 2 

ISR-jet 

Pt (ISR-jet) > 550 GeV 

tops 

m(top-candidates) < 250 GeV 

alignment 

> 0.55, |A(/)(ISR-jet,|;T)| > 2.95 

f^p/\/H t 

optimized window (minimally fix /\/ Hp > 3 GeV* A) 


TABLE I: Summary of reconstruction cuts. 


level identity (fe-j^t or light-flavor jet) based on the above efficiencies. Mistags of light-flavor 
jets are not incorporated, nor are backgrounds with less than two heavy-flavor partons 
in the hard event (see Appendix 0. Light-flavor mistags are of subleading importance 
for both the stop signal and top backgrounds. For W/Z+]ets and especially multijets, a 
complete analysis with light-flavor mistags requires extensive simulation, which we have not 
undertaken. However, we do not expect this omission to have significant impact on the 
validity of our background estimates. As a specific corroborating example, we refer to the 
detailed background composition of the Higgs search (W/Z)H —)■ {W/Z){bb) |36], in which 
the W/Z+]Qts backgrounds are dominated by events with two truth 6-jets."^ 

The vector is modeled in two ways, with the hopes of bracketing true performance. 
The first model is truth which though simple to define in simulation is clearly overly 
optimistic. The second model is firp = — 'YhjPTij)- This definition is technically somewhat 
pessimistic because it does not account for corrections that may come from activity not 
clustered into jets. However, we re-emphasize that pileup effects have been incorporated 
into the jet energy resolutions. For both definitions, is not allowed to point along 

the pt of any of the leading three jets, with a requirement |A0| > 0.55. In practice, such 
a cut is used experimentally to avoid fake from under-measured jets, as well as real 
from heavy flavor decays inside of jets. Within our own multijets samples, the cut is still 
somewhat advantageous when using The advantage with truth ^p is minor, but we 
continue to apply the cut to maintain consistency and a higher degree of realism. 

Identification of the ISR jet exploits the kinematics of top decay in a simple way. For a b- 
quark produced in a hadronic top decay, adding in either of the quarks produced in the sister 
LF’s decay will produce a subsystem with a mass less than m*, and more specifically less 

^ To give some rough sense of accounting, the “penalty” for QCD to produce a pair of hard, well-separated 
heavy quarks from a gluon splitting is 0(as/7r), which is overall percent-scale. This easily beats the chances 
of a double-mistag of truth light-flavor jets, which is 0(10“'*). For single-mistag events containing one 
6-quark at the hard event level, the 0(10“^) mistag would need to be combined with the very small 6 
PDFs. (Practically such events are paying both Os/tt a,nd the mistag rate.) 
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FIG. 1: Kinematic distributions of stacked backgrounds and some example signal points for the 
.^T-based analysis, with cumulative cuts. The baseline cuts include the lepton veto, jet counting, 
ISR-jet pt cut, and a cut fix/\/H t > 3 GeV^/^ used to define the simulation samples. 


than TTif — ~ 153 GeV at leading-order with narrow W. These ineqnalities continue 

to hold even when the top is below its mass-shell, as the kinematic boundary only becomes 
lower. The leading two &-jets in the event are taken to be the 6-quark candidates. A list 
of remaining jets in the event is formed which satisfy m{b + j) > 200 GeV for both 6-quark 
candidates. The highest-pT jet from this list is then the ISR candidate. Only events with 
Pt (ISR-jet) > 550 GeV are kept in our analysis. 

Individual top quarks are reconstructed using a procedure borrowed from [5] . Excluding 
the two leading 6-jets and the ISR-jet candidate, the two closest jets in the p-0 plane are 
added to form a “hT boson.” This in turn is added to the closest 6-jet to form a “top 
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FIG. 2: Distribution of .^^^-significance for stacked backgrounds and some example signal points, 
for the f^T analysis (left) and the truth-analysis (right). All other cuts have been applied. 


quark.” The procedure is then repeated amongst the remaining jets and 6-jet. In the 
absence of smearings and combinatoric confusions, both top-candidates constructed in this 
manner would satisfy m ~ if on-shell, and m < if off-shell. We make a somewhat 
looser demand of m < 250 GeV. The main purpose of this cut is to reduce multijet and 
lT/Z-|-jet backgrounds, which tend to reconstruct higher masses with a very broad tail. 

Finally, we employ the relation in Eq. which, as per |30], we decompose into angle and 
magnitude. For the angular component, a strong anticorrelation between the ISR-jet and 
directions is demanded: |A0(ISR-jet,.^T)| > 2.95. For the magnitude, we expect that the 
signal will be approximately equal to (ISR-jet) x (m^/mj). Because of the interplay of 
the hard (ISR-jet) cut and the rapidly-falling production pt distributions, the signal will 
appear as a localized bump in Raw can serve as an adequate discriminating variable 
here, but we instead use the “signihcance” ratio ^t/\/Ht, which we hnd to be slightly more 
efficient at separating signal from background for the lower-mass signals. 

Table H] summarizes the complete set of cuts. Figs. [T] and show distributions of several 
of the discriminating variables for backgrounds and some example signal points, illustrating 
the cumulative purihcation of the signal. Note that, to maintain efficient Monte Carlo 
generation, a cut of fix! a/ Ht > 3 GeV^/^ has been applied to dehne a baseline reconstructed 
sample (regardless of the hnal dehnition used). 

The analysis thus dehned, we scan through the model space of the stop-neutralino mass 
plane, with hner-grained steps near the top compression line (of order the top quark width). 
The hnal / a/ Ht window is optimized per sample to maximize the naive statistical signih¬ 
cance S/\/S + B. We dehne exclusion threshold as 5/ -s/S + B = 2, and discovery threshold 
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FIG. 3: Existing Run 1 limits from ATLAS and CMS, and projected 300 fb“^ discovery and 
exclusion sensitivities for our analysis. The truth-analysis (not shown) yields very similar 
exclusion contours, but somewhat stronger discovery contours at lower masses. Note that our 
simulation grid does not extend all the way down to the W compression line ~ + mw 

nor below, where the decay kinematics transitions to four-body. (We also do not indicate existing 
exclusions in that region. For the stealth exclusions, see Fig. [^) 


as S/y/B = 5. Fig. shows our final exclusion and discovery contours for 300 fb indicat¬ 
ing a near complete closure of the current compression line gap. Fig. shows the luminosity 
required to achieve exclusion-level sensitivity along the compression line. While our simu¬ 
lations are done under Run 2 & 3 conditions, we have also naively extrapolated as far as 
the HL-LHC luminosity of 3 ab“^. We include as well in Fig. a scan of the signal and 
background rates at 300 fb“^ along the top compression line. This indicates S/B ^ 1 over 
most of the range that we study, suggesting good resilience to systematic errors, which we 
have not attempted to estimate. Finally, in Fig. |^we provide a closer view of the exclusion 
sensitivity near the stealth point, via a series of scans over at hxed neutralino masses. 


III. DISCUSSION AND OUTLOOK 

Fig. suggests that our proposed search strategy can access stops along the top com¬ 
pression line beyond 400 GeV at discovery-level signihcance, and perhaps up to 550 GeV at 
















Lint needed for exclusion along top compression line 



FIG. 4: Luminosity required for exclusion sensitivity along the top compression line for both our 
truth-(dashed) and fix (solid) analyses, assuming 13 TeV and Run 2 & 3 pileup and detector 
conditions. (Projections beyond 300 fb“^ are naive extrapolations, not using HL-LHC conditions.) 



FIG. 5: Signal (red) and background (purple) counts along the top compression line for 300 fb ^ 
for both our truth-(dashed) and fix (solid) analyses. 


exclusion-level significance, over the current phase of LHC running. These numbers already 
start to approach what was done for non-compressed stops at Run 1. However, unlike those 
searches, for us the sensitivity is maximized on the top compression line. This complemen¬ 
tarity is made possible by focusing on the unique kinematic conhgurations that start to open 
up at Run 2. It is rather remarkable that the sensitivity gap at the top compression line, 
which has become a modern benchmark of difficulty in new physics searches, can be shrunk 
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FIG. 6 : Projected 300 fb“^ exclusion sensitivity around the stealth point, scanning for sev¬ 
eral fixed <C mt, for our fix analysis. The gray shaded region indicates the excluded stop 
masses from the dedicated ATLAS search for stealth stops [T] . The truth-.^ 7 ’ analysis (not shown) 
yields very similar significances in the dip near the top compression line, but up to 0 ( 1 ) higher 
significances away from it. 


so dramatically by simply going to a higher collider energy. Fig indicates that the gap 
will start to close already with a data set comparable in size to Run 1, which should be 
achievable before the end of 2016. 

On the low side, our search very closely surrounds the stealth point (mp m^) = {rrit, 0), as 
indicated in detail in Fig. In fact, we have found that the exclusion-level contour there de¬ 
pends only moderately on whether we use truth or fJx, though Fig. [^illustrates that this 
choice does strongly effect the S/B there. We emphasize the caveat that we have not folded 
in systematic errors. Ultimately, the major question is how well the multijet background 
can be controlled and modeled. Given this uncertainty, it is difficult for us to make very 
concrete statements near the stealth point. But following the discussion in the introduction, 
it seems highly likely that multiple search strategies will come into play. Even the present 
state-of-the art searches based on tt cross section and spin correlation measurements Pi 
already overlap with our projections, completing the coverage at exclusion-level. 

Our search is also very effective at covering large portions of the three-body region. 
While our simulated model points do not extend below the W compression line at ~ 
+ nib + fnw, and into the four-body region, it seems quite likely that we even continue 
to have some coverage there. This leaves open the possibility of linking up with monojet 
and other searches in that region. (See as well [2S] for a recast of a soft dilepton search at 
7 TeV that already makes some surprising inroads there.) An approach that requires fewer 
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jets and looser hadronic top reconstructions would also likely be fruitful, a possibility that 
we save for future work. 

More generally, we have only very coarsely optimized our analysis, hrst by hxing most 
of our selection criteria by-eye on a small subset of model points, and then by selectively 
scanning over only our hnal ^t/\/Ht window. With the principle proven, a more carefully 
optimized suite of cuts would certainly achieve better results, especially for the stealthier 
model points. Breaking the search into more analysis regions, e.g. binned over p 7 ’(ISR-jet) 
(or £t over multiple variables), could also be beneficial. 

An obvious further extension of the analysis includes HL-LHC, with up to 3 ab“^ of 
luminosity. The very high pileup would likely be a major concern there, as the rate of fake 
jets rises signihcantly, and the resolution on further degrades. Certainly, pushing further 
into the stealth region will be difficult, although the much higher event rates may allow for 
more highly-crafted cuts. On the high-mass side, if we naively extrapolate up our 300 fb“^ 
analysis as per Fig. we hnd discovery (exclusion) reach extending to about 800 GeV. 
Along similar lines, projections for a 100 TeV proton collider are also interesting to pursue. 
However, as we ultimately scan up to mi 3> mj, we effectively return to the fully compressed 
situation m^. ~ m^. All of the compression lines may then practically blur together using 
more standard “monojet”-|-.^T style searches, perhaps supplemented by the additional “soft” 
activity from the decays. Such an analysis has been carried out in EH, hnding sensitivity 
to compressed stops up to multiple TeV using the dilepton channel.® Finally, all of our 
results readily generalize to those classes of fermionic top-partner models that exhibit either 
a conserved or approximately-conserved parity, and contain a neutral “LSP” boson which 
plays a role kinematically identical to [38]. The only major difference relative to stops, 
from the perspective of our analysis, is their approximately six times larger cross section at 
a given mass, yielding commensurately stronger sensitivity. 

In conclusion, natural supersymmetry poses some interesting phenomenological chal¬ 
lenges, as evidenced by the enduring gaps in coverage of one its simplest incarnations: 
an NLSP stop and LSP neutralino. While limits continue to push upward in mass in the 
favorable parameter regions that readily provide lots of ^t, we have seen here that an ap¬ 
propriately constructed analysis at the upgraded LHC, along the lines suggested in [30], can 
qualitatively extend sensitivity to this model into the more difficult compressed regions at 
lower masses. Combined, these approaches will leave very little “natural” parameter space 
unexplored. With its next major phase in progress, the LHC appears poised to provide us 


® If we naively scale the energies and cross sections from the existing monojet-I-searches for fully 
compressed stops ms] from an 8 TeV machine to a 100 TeV machine (without running the PDFs), 
we would expect an exclusion of (260 GeV) x (100/8) ~ 3 TeV after accumulating a luminosity of 
20 fb“^ X (100/8)^ ~ 3 ab“^. Suggestively, this coarse estimate is very close to that of [37] on the top 
compression line. 
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with a much more comprehensive perspective on the possible role of supersymmetric top 
quarks in Nature. 


Note added: While this paper was nearing completion, [SH] appeared, which has signih- 
cant overlap with our results. Their proposed Rm variable (a very close variant of what was 
originally proposed in [30]) is highly correlated with the ^t/\/Ht variable that we use here, 
and in general with any variable proportional to in the presence of a hard ISR-jet pt cut. 
There are a number of other differences in our analysis strategy, which lead to a higher S/B 
with comparable formal statistical significance, and somewhat different sensitivity contours. 
We also pay additional attention to the approach to the stealth region and the possible role 
of resolution. However, we do not make a dedicated study around the W compression 
line. 
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Appendix A: Event Generation 

Our event generation is performed using MadGraph5_aMC@NL0 jlO] at 13 TeV and showered 
with PYTHIA 6 mi, using leading-order matrix elements (without K-factors). We set the 
top quark mass to 173 GeV, and width to 1.5 GeV. 

For our signal samples, we choose mostly-right-handed stop and mostly-Bino neutralino 
(spin effects on our all-hadronic analysis are expected to be modest). Most samples are 
generated as tt*j, with only a parton-level cut of 400 GeV on the accompanying jet. Both 
stops are decayed using three-body phase space i —>■ Wbx^, regardless of mass point, which 
is crucial for modeling the kinematic transition at the top compression line. A complete 
decay chain is therefore, e.g., tl > W+ b nl, W+ > j j. The stop width for each model 
point is computed separately using 1 —)■ 3 parton-level decay simulations. A subset of models 
along the compression line have been simulated over their full production phase space, using 
fc-p-MEM matching with a threshold of 100 GeV. Perhaps unsurprisingly, the events passing 
our final selections are highly dominated by the Ij subsample, and are in close agreement 
with our simple unmatched simulations. Similarly, we hnd very low relative pass rates for 
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decay modes other than all-hadronic. 

The backgronnds are generated as follows.® Onr tt sample is matched np to one (two) 
jets for all-hadronic (partially leptonic or r) decays, again nsing a 100 GeV threshold. We 
also generate tiW and ttZ matched np to one jet. For W/Z+]ets and mnltijet backgronnds, 
we concentrate on prodnction with at least two heavy qnarks (bottom or charm) in the hard 
event. Becanse of the difficnlties of compnting very high-mnltiplicity matrix elements, we 
mainly use the parton shower to generate extra partons, and do not employ any matching. 
The W/Z+]ets sample specifically starts with W/Z (decaying to Iv, rz/, or vv) plus three 
hard partons, while the multijet sample starts with four hard partons. We have also cross¬ 
checked the multijets against AlpGen |l2] samples, generated with identical criteria. For each 
sample we impose cuts at the parton level that treat the b and j partons democratically, 
requiring prij) > 15, AR{j,j) > 0.4 (where j here includes b) as well as a pt cut on the 
hardest jet of 350 GeV. 
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